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Abstract 


Marginal short-period contact binaries are important to understand as they pose a different physical scenario than 
the predicted theoretical model based on the thermal relaxation oscillation mechanism due to their shallow degree 
of contact. Here we present the optical and X-ray studies of a contact binary source RW Dor using the Transiting 
Exoplanet Survey Satellite (TESS) and XMM-Newton telescopes. For the first time we report the varying 
O’Connell effect and explain the asymmetry with a spot model. Based on the new times of minima, we make a 
robust estimate of the orbital period of the third body at 47.01 ~ 0.52 yr with an eccentricity e = 0.21. We show 
that the period-decreasing trend observed in O — C variation can be explained by both conservative mass transfer 
from primary to secondary and AML via stellar wind. The X-ray luminosity exhibited by RW Dor did not vary 
significantly on three different occasions and was found to be about 3.34 x 102? ergs '. Assuming that the 
quiescent X-ray emission is emitted from an undisturbed loop structure, the loop size is estimated to be 
0.6-1 x 10? cm which is <Alfvén radius ra ^ 8 x 10'? cm. 
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1. Introduction 


Contact binaries are a subclass of close binary systems 
whose primary or secondary components have filled their 
Roche lobes and share a common convective envelope (CCE) 
with a varying degree of fill-out factor (Kopal 1959; 
Lucy 1968). The fundamental parameters of these binaries 
provide us with a better understanding of some of the most 
interesting physical processes, such as mass transfer, 
magnetic cycles, strong interactions between stars and stellar 
mergers (see Qian et al. 2017 for a detailed review). There 
exist two types of contact binaries, i.e., A- and W-subtypes 
based on the temperatures and masses (Binnendijk 1970). In 
the W-subtype, the primary eclipse is due to the more 
massive cooler component transiting the less-massive hotter 
component, and it is the opposite scenario for the A-subtype. 
Based on the long-term periodic studies, it also draws a 
picture of the presence of a third body in many such systems 
which can explain the low angular momentum configuration 
often observed in these systems. The optical imaging by 
Ruciński (1969) clearly indicates the presence of a third body 
in most contact binaries, explaining their pivotal dynamical 
role in the evolution of these systems. This was further 
confirmed by the O — C diagram-based studies carried out on 
those binary systems. 

Contact binaries are known to exhibit the O'Connell effect 
(O'Connell 1951) which is the asymmetry observed in the 
photometric light curve generally attributed to magnetic activity 


over the stellar surface. The synchronous fast rotation of 
common envelopes is somehow associated with the origin of 
X-ray emission (e.g., Gondoin 2004a). These systems exhibit 
high chromospheric activity along with high coronal emission 
and thus are strong X-ray emitters (Carroll et al. 1980; 
Huenemoerder et al. 2006; Hu et al. 2016). It has been observed 
that the strength of X-rays is correlated to the orbital period and 
also the temperature of binary companions (Stepien et al. 2001). 
It was noted that the massive component is responsible for the 
X-ray emission through their magnetic activity, e.g., VW Cep 
and Y Y Eri (Vilhu & Maceroni 2007). The Advanced Satellite 
for Cosmology and Astrophysics (ASCA) spectra displayed the 
characteristics of a corona and found two-component model 
temperatures of 7 x 106 К and 22 x 10$ K along with a flux of 
1 x 10 !'ergem 287! and VW Cep appears to show polar 
spots based on Doppler imaging maps (Hendry & 
Mochnacki 2000). It has been found that magnetic loops in 
rapidly rotating low-mass stars would sweep to the poles 
(Buzasi 1997). Moreover, contact binaries display a super- 
saturation effect (McGale et al. 1996; Chen et al. 2006) and it 
has been argued that this effect is due to the large loops 
becoming unstable by the Coriolis force and being extended 
along the poles (Jardine & Unruh 1999). In this model, the loops 
are large and extended but Stepien et al. (2001) showed that the 
loops could be compact when compared to the stellar radius and 
spatially associated with the equatorial region of the stellar 
surface. The study of O — C curves of overcontact binaries by 
Tran et al. (2013) using the Kepler satellite gave anticorrelated 
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O — C curves for the primary and secondary eclipses that they 
interpreted as the outcome of the longitudinal movement of 
starspots on the stellar surface. A detailed study of O — C curves 
indicates that the movement of a stellar spot is associated with 
the rotational motion of the individual star, and is asynchronous 
in 5096 of the cases. 

RW Doradus (RW Dor) is an interesting binary system with 
a short orbital period of 0.2854 day exhibiting G4/5 spectral 
type. The light curve displays a difference in the depth of 
minima and shows asymmetry in the light curve which was 
modeled as a hot spot in the massive component around the 
joining of the two stars (Marton et al. 1989). However, Kaluzny 
& Caillault (1989) did not find any significant asymmetry in the 
light curve, rejecting the idea of a spot over either of the stellar 
components. Marino et al. (2007) did observe asymmetry in the 
light curve but a strong O'Connell effect can be ruled out due 
to the large scatter in the maximum phase of the light curve. 
They also observed an orbital period that was decreasing with a 
rate of AP/P = —6.3 x 10 !'. In addition, Deb & Singh 
(2011) did not find any asymmetry in the light curve which was 
confirmed by the absence of a third light. Sarotsakulchai et al. 
(2019) did not find any observational signature of asymmetry 
in the light curve or O'Connell effect, strongly indicating no or 
a low level of stellar activity. They also found weak evidence 
of a decreasing period along with a cyclic trend in the O — C 
curve, suggesting the possible presence of a third body. The 
radial velocity measurements found a spectroscopic mass ratio 
р = 0.68 with a similar spectral type exhibiting features of a 
КІ type star (Hilditch et al. 1992). Later, a similar qsp = 0.63 
was reported by Duerbeck & Rucinski (2007) who found 
V, — 25 kms. Overall, these studies suggest that no strong 
activity has been observed in RW Dor with weak evidence of 
periodic variations in the respective light curve studies. 

In the present work, we focus on the modeling of the light 
curve observed by the Transiting Exoplanet Survey Satellite 
(TESS) which displayed an O'Connell effect for the first time. 
Detailed studies have been performed to find a robust estimate 
of the periodic variation using previously published times of 
minima along with new ones with TESS and the All Sky 
Automated Survey (ASAS). We also studied the X-ray 
emission on different occasions observed by the XMM-Newton 
observatory to check the variations in the X-ray flux in the 
observed duration. 


2. Observation and Data Analysis 


TESS was placed at a stable orbit in a 2:1 resonance with the 
Moon following an elliptical orbit of 13.7 days (Ricker et al. 
2015). The goal of the satellite is to look for exoplanets orbiting 
bright stars using the transit method. TESS observes the 
selected fields with a 10.5cm optical telescope in the 
wavelength range of 600-1000 nm at 2100 pixel '. Its field 
of view is 2 x 96 deg? and scans the sky for 27 days. RW Dor 
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was observed by TESS on two occasions, i.e., 2018 October 19 
and 2018 December 15, with a cadence of 2 minutes. Pre- 
search Data Conditioning Simple Aperture Photometry 
(PDCSAP) data were used to study the light curve variation, 
and data were downloaded from the Mikulski Archive for 
Space Telescopes (MAST).! 

We have also studied the archival data of the XMM-Newton 
observatory taken on three different occasions for the source 
RW Dor (a = 5" 18" 3275, 5 = —68° 13/ 327) that was first 
observed on 2009 May 16 (ObsId.:0602980201) for 20.9 ks for 
which only data from the EPIC MOS1 camera were available. 
The second observation was on 2010 December 11 
(ObsId.:0650020101) for a duration of 24.1 ks and was 
observed by both the EPIC MOS(1&2) and PN cameras. The 
third observation was on 2013 October 22 (Obs. 
ID:0723650201) for a duration of 20ks and only EPIC 
MOS1&2 data were available. We have selected a circular 
source region and a background region having radius 50" for 
our first and second observations, while for the third 
observation we have used a circular region with radius 35" 
for both source and background regions. The source and 
background spectra were obtained using the expression 
"PATTERN « 12" for MOS data and "PATTERN « 4" for 
PN data using SAS version 19.1.0. We have obtained the 
redistribution matrix file and the ancillary response file utilizing 
rmfgen and arfgen tasks respectively. The spectral analysis was 
done employing XSPEC v 12.12.0 (Arnaud 1996). 


3. TESS Detection of O'Connell Effect and Light 
Curve Solution 


Previous optical studies of RW Dor have not reported any 
strong signature of stellar activity, i.e., O'Connell effect (for 
more details see Sarotsakulchai et al. 2019). The TESS 
PDCSAP light curve clearly exhibits the variable O'Connell 
effect (Figure 1, top panel). Figure 1 is obtained by taking the 
difference between the consecutive maxima in the observed 
light curve. During the first observation, which lasted for 
25 days, a weak signature of O'Connell effect was observed, 
however, in the second span of observation, the effect was 
relatively evident. It can be seen that initially the O'Connell 
effect was not observed and gradually it started appearing. We 
fitted the phased light curve with and without the O'Connell 
effect (Figure 1, bottom panel). The light curve solutions were 
obtained using the Wilson-Devinney code-based program 
PHOEBE (Prsa & Zwitter 2005). 

The following methodology was invoked to get the best 
solution for the light curves (for more details see Sriram et al. 
2016, 2017, 2018). The primary component temperature was 
fixed at 5560 K based on a recent study by Sarotsakulchai et al. 
(2019) who constrained this value based on B — V — 0.69. We 
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Table 1 
Photometric Solutions Obtained for RW Dor Using TESS Data 

No Spot With Spot 
Parameter 
A, = А 0.50 0.50 
81782 0.32 0.32 
T\(K) 5560 5560 
T(K) 5200 + 69 5200 + 16 
q 0.6147 + 0.0063 0.6147 + 0.0026 
i^ 76.8334 + 0.03 76.3580 + 0.12 
Q=% 3.0583 + 0.0055 3.0816 + 0.0034 
rı (pole) 0.4042 0.4005 
rı (side) 0.4294 0.4247 
rı (back) 0.4626 0.4563 
r(pole) 0.3195 0.3195 
rj (side ) 0.3343 0.3343 
r2 (back) 0.3672 0.36798 
Lı /(Li+L2) 0.4459 0.4362 
Spot colatitude (0?) Ue 3243 
Spot longitude (c?) EE 40.3 + 2.1 
Spot radius (r^) er 15.2 + 3.6 
Temp. ratio Ut 0.69 + 0.08 
Xw(o — с)? 0.0061 0.0052 


fixed the gravity darkening coefficient gı = g2 = 0.32 
(Lucy 1967) and adopted albedos A; = А5 = 0.5 for both 
components (Ruciński 1969). Bolometric limb darkening 
coefficients (logarithmic law) were adopted from the table 
determined by van Hamme (1993). A circular and synchronous 
(F=1) orbit was assumed. Initially, four parameters were 
adjusted, 1.е., the temperature of the secondary component (75), 
orbital inclination (7), the dimensionless potentials of the 
primary and secondary components ((); = Qı), and the 
luminosity of the primary star (L4). The mass ratio was fixed 
at 0.63 based on the spectroscopic observations (Duerbeck & 
Rucinski 2007) and later freed along with other parameters to 
get the best solution. The best-fit solution resulted in a mass 
ratio а = 0.614 + 0.006 and an inclination i — 77? + 0.12. The 
best-fit results are shown in Table 1 and the obtained values are 
similar to those reported in previous studies (Sarotsakulchai 
et al. 2019). In Figure 1, the middle panel displays the best fit 
using the black line. It is clear from the fit that a no-spot 
solution produces more residuals around phase 0.75. Later, we 
adopted a cool spot over the primary component and continued 
the fitting. The cool spot model has four parameters: the 
colatitude (0), longitude (©), and radius of the spot (r,,4) and 
spot temperature ratio (Tspot/ Tstar). The best-fit parameters аге 
given in Table 1 after including a cool spot and finding the size 
of the stellar spot to be 15° at a colatitude of ~32° during the 
Heliocentric Julian Date (HJD) 2458489.2509 to HJD 
2458489.3079 and the temperature ratio is ~0.65. In order to 
produce a robust the solution, we varied the primary 
component’s temperature from 5400 to 5800 K and performed 
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the fitting. The lowest residuals were noted to be the best 
solution which is displayed in Table 1 (Figure 1, bottom panel). 


4. Period Variation Studies 


Previous studies by Sarotsakulchai et al. (2019) observed a 
period-decreasing trend along with weak evidence of sinusoidal 
variation, indicating the possible presence of a third body with 
an orbital period of ~50 yr. Many times of minima were older 
than 1940 and hence we used TESS data to make the study of 
the О — C variation more robust. The new times of minima 
along with others are logged in Table 2. The TESS time units are 
converted to HJD using the utility mentioned in the link https:/ / 
astroutils.astronomy.osu.edu/time/index.html. In order to com- 
pletely determine the various parameters of the third body 
assuming cyclic variation, we applied Zasche's code (Zasche 
et al. 2009) in MATLAB which implements the Simplex method 
in order to arrive at best-fit solutions. The following equation 
was fitted to the overall O — C diagram in Figure 2 using a red 
line and a blue line displays only the period-decreasing quadratic 
term. 


Min. = JD, + PE + QE + 227707 
l-e? . | 
x | ——— — — sin(v + шз) + ез sin(ws) |, (1) 
1 + escosv 


where a2 sin i, c, ез, v and шз are the projected semimajor axis, 
speed of light, eccentricity, true anomaly of the binary orbit 
around the center of mass of the triple system and the longitude 
of the periastron, respectively. Table 3 displays the best-fit 
values for various parameters for different times of minima in 
the initial three columns. 

To make the analysis of the third body more robust, we 
varied the eccentricity of the third body, i.e., ез = 0-1 and 
allowed other parameters to vary, and noted the lowest residual 
at ез = 0.21. Later, we kept ез = 0.21 which was freed 
subsequently and obtained the best solutions (see Table 3). 
To confirm the third body orbital period, we took different 
times of minima as shown in Figure 2 and we found that the 
orbital period varied from 46 to 54yr. The best period 
P3 = 46.90 + 0.73 was obtained for JD, = 2458435.0148 
based on the sum of the residuals. We also performed a similar 
analysis by considering the times of minima observed after the 
year 1980. Figure 3 shows the best fit and we noted a period of 
Рз ~ 26.55 + 1.13 yr (see Table 3 last column). 

To make the result even more robust, we randomly selected 
epochs after the year 1985 (see Table 2) and estimated the third 
body orbital period. The distribution of the period is depicted in 
the histogram (Figure 4, top panel) and it is clear that the peak is 
observed around 47 yr. A Gaussian profile is fitted to this peak 
which is displayed with a dashed line in Figure 4 (top panel) and 
the best fit resulted in a mean period of Р; џ= 47.01 and 
о = 0.52 yr. The higher periods were sparse and resulted due to 
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Figure 1. (a): Top panel displays the flux difference around the consecutive maximum and the dashed line at zero corresponds to the no O'Connell effect. It is clear 
that the O'Connell effect varies during the TESS observations. Time is in units of Barycentric Julian Date (BJD). (b): Middle panel displays fitted photometric 
solutions to the phased light curve along with a no spot solution (black line) and a spot solution (red line). (c): Bottom panel depicts the residual of the fits for different 
primary component temperature (T1) (see text for details). 
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Table 2 
All Available Times of Minima for RW Dor 
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(O-C) 


—0.229449 
—0.23324 
—0.045714 
—0.008073 
—0.007215 
—0.005538 
—0.006511 
—0.005094 
—0.004547 
—0.00613 
—0.002273 
—0.002396 
—0.002512 
—0.00262 
—0.00325 
—0.003938 
—0.004337 
—0.000516 
— 0.001663 
—0.005429 
— 0.003588 
—0.006553 
— 0.005283 
— 0.004865 
—0.004643 
—0.004136 


—0.002234 
—0.002423 
—0.002812 
—0.002201 
—0.00279 
—0.002079 
— 0.002068 
—0.002157 
—0.002376 
— 0.002465 
—0.002258 
—0.002447 
—0.002562 
—0.002451 
—0.00244 
—0.002629 
—0.002518 
—0.002607 
—0.002796 
—0.002785 
—0.002374 
— 0.002963 
0.000648 


(0 — С) 


0.00526872 
—0.0164409 
0.0357531 
0.00296067 
0.00373861 
0.00480244 
0.00369319 
0.00486953 
0.00541339 
0.00359229 
0.00659192 
0.00645324 
0.00321103 
0.00307041 
—0.00346362 
—0.00483576 
—0.00541813 
—0.00177321 
—0.00291681 
—0.00667067 
—0.00373681 
—0.00536357 
—0.00363601 
—0.00290571 
—0.0025953 
—0.00192461 


0.000939127 
0.000751621 
0.000364116 
0.00097661 
0.000389105 
0.0011016 
0.00111409 
0.00102659 
0.000814073 
0.00072657 
0.00093906 
0.000751558 
0.000639055 
0.000751553 
0.000764052 
0.000576551 
0.00068905 
0.00060155 
0.00041405 
0.00042655 
0.00083905 
0.00025155 
0.00386405 
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Figure 2. (a): The top panel displays the best-fit O — C variations using LITE 
represents the LITE solution. (a^: The corresponding right panels show only the 


Table 2). 


the older epochs. A very similar analysis was performed for the 
times of minima after the year 1980 which are relatively more 
robust. Figure 4 (bottom panel) displays the histogram along with 
a Gaussian fit with a dashed line. The best fit resulted in a mean 


period of Р; u — 27.22 and с = 3.12 yr. 


along with the residuals. The blue dashed line is a parabolic fit and the dashed line 
LITE fit to the data. Other panels display the same with different times of minima (see 


5. X-Ray Spectral Analysis 


We explored whether the X-ray emission of RW Dor varied 
during a span of four years as observed by the XMM-Newton 
telescope. It was done because previous studies did not observe 
any stellar activity in the form of the O'Connell effect 
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Figure 3. (a): Top panel: O — C variations along with the LITE solution (red line) and quadratic fit (blue line). (b): The bottom panel displays the cyclic variation. 


Table 3 


Third-body Solution Based on O — C Diagram for Different Epochs 


Parameter 


A, semi. amplitude (days) 

Рз, period (yr) 

ез, eccentricity 

шз, longitude of periastron passage 
T,, time of periastron passage 
алә Sin i, (au) 

Q 

(M3), (Mc) 

M3(Mo)i=90° 

M3(Mo)i=60° 

Ms (Mo)i=30° 

Sum of the square residuals 


HJD = 2457112.6116 


0.0099 + 0.0013 
54.51 + 3.55 
0.21 + 0.49 
0.00 + 78.11 
2440313.946 + 5118.252 
1.75 + 0.23 
—0.0998 + 0.00002 
0.0018 + 0.000016 
0.160 + 0.00051 
0.187 + 0.00060 
0.346 + 0.0011 
0.0116 


HJD — 2454037.5995 


0.0045 + 0.0018 
40.75 + 2.88 
0.21 + 0.71 
9.183 + 117.88 
2443798.757 + 4330.103 
0.801 + 0.328 
-0.0938 + 0.00001 
0.000311 + 0.000001 


0.085 + 0.00014 

0.099 + 0.00017 

0.179 + 0.00032 
0.0135 


Note. The last column is obtained only for times of minima after the year 1985 (see text). 


HJD = 2458435.0148 


0.0049 + 0.0003 
46.90 + 0.73 
0.21 + 0.12 
0.81 + 22.42 
2442530.250 + 1058.911 
0.8751 + 0.057 
—0.111 + 0.0001 
0.00030 + 0.000001 
0.0852 + 0.002 
0.100 + 0.001 
0.177 + 0.002 
0.0016 


HJD = 2458436.4387 


0.0045 + 0.0004 
26.55 + 1.13 
0.22 + 0.16 
211.60 + 27.57 
2454034.424 + 735.439 
0.7953 + 0.0740 
—0.218 + 0.0005 
0.00071 + 0.0000040 
0.1147 + 0.0002 
0.1336 + 0.0003 
0.2427 + 0.0005 
0.0075 
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Figure 4. (a): Top panel: Histogram of the third body period obtained from the LITE solution with different epochs after the year 1985. The dashed line is a best-fit 
Gaussian profile with и = 47.01 апас = 0.52. (b): Bottom panel: Same as above but with the LITE solution considering the data after 1985. The dashed line is a best- 


fit Gaussian profile with и = 27.22 and о = 3.12. 


(Sarotsakulchai et al. 2019). We fitted the X-ray spectra of 
MOS1, MOS2 and PN as shown in Figures 5—7. First, we tried 
with a simple power-law model to fit the spectra but the 
observed x were noted to be 22. Then we unfolded the 
spectra with wabs*bremsstrahlung and wabs*bbody models. 


The hydrogen column density was fixed to 5 x 10 cm ? 
based on Dicky & Lockman HI in the Galaxy (Dickey & 
Lockman 1990). In the case of the bremsstrahlung model, 
it was found that electron temperature was around 
kT. = 0.60keV with a poor fit (Table 5). The blackbody 
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Figure 5. The top panel shows the unfolded spectrum along with the fit (thick line) for EPIC MOS1. The bottom panel displays the residuals. The reference number is 


ObsID.0602980201 and observation date is 2009 May 16. 


temperature was about KT ~ 0.17 + 0.02 on all the occasions 
(Table 6), not displaying any variation. 

Later we attempted a plasma model, i.e., apec model (Smith 
et al. 2001), which has two main parameters, viz. plasma 
temperature kT, and abundance Z (Zo). Initially, a single apec 
model resulted in a high x?, and hence a two-component apec 
model was used. The results are shown in Table 7 for Mos 1 
and 2 and PN for different years. The kT, = 0.24 + 0.08 keV 
and kT, = 0.93 + 0.15 keV for the year 2010 along with flux 
1.61 x 10 P ergcm ?s^!. In general, if the source would 
have been active then the active coronal X-ray emission would 
have three or four components which was not evident from the 
present studies. The observed coronal temperatures strongly 
indicate the presence of a quiescent corona in RW Dor and 
found no significant variation from 2009 to 2013. The X-ray 
luminosity was estimated based on the distance determined 
from the Gaia parallax, i.e., d — 123.57 pc (Table 7). 


6. Results and Discussion 
6.1. Conservative Mass Transfer 


For the first time, we report clear evidence of the variable 
O'Connell effect in RW Dor. Based on the photometric 
solutions, we observed that a cool spot of size 15? at a latitude 
of 60? on the primary component is necessary for the best fit. 
The photometric solution mass ratio was found to be g — 0.61 
with an inclination of 77? which is well in agreement with 
previous studies (Sarotsakulchai et al. 2019). Based on 
previous times of minima along with TESS and ASAS data, 
we found that the orbital period of RW Dor is decreasing at a 
rate of —2.80 x 10 ? day yr! along with a cyclic variation. 
This suggests that material is transferring from the primary to 
the secondary component, driving the mass ratio of the system 
to increase. Assuming a conservative mass transfer rate from 
primary to secondary, it can be estimated from the following 


equation 


m 3M, | : | | (2) 
Р M M» 

Assuming M, = 0.52 Mo and М = 0.82 Mo, M, was found to 
be 4.64 x 1078 M. yr !. The timescale of mass transfer is 
M,/M,= 11.2 Myr. The mechanism of angular momentum loss 
(AML) due to magnetic stellar wind has a tendency to decrease 
the orbital period during a timescale of 8.6 Myr (Sarotsakulchai 
et al. 2019). Based on the results of the present study, the time 
P/ P=10 Myr, which is close to the value obtained from the 
AML mechanism. This indicates the magnetic stellar wind can 
also drive the decrease in the orbital period as the spot activity 
has been initiated in the source. These systems often go through 
thermal relaxation oscillation (TRO) along with the AML 
mechanism, whose period is of the order of a few million years 
(Qian 2001). As the period continues to decrease, the degree of 
contact between the components must increase. In the present 
study, we observed evidence of a cool spot on the primary 
component, suggesting that variable activity has been initiated 
which was not observed in earlier optical photometric or 
spectroscopic studies. This stellar activity would have triggered 
the AML and contributed to the decrease in the orbital period in 
RW Dor. In general, as the degree of contact increases it 
reduces the AML as the convective envelope probably veils the 
activity (e.g., Qian 2001). However, in the present scenario, it 
is not the case because of the appearance of the activity and this 
system has a shallow degree of contact (f~ 10%). Future 
observations are necessary in order to confirm the results as 
well as to track the stellar activity of the source. 


6.2. Non-conservative Mass Loss in RW Dor 


Justham et al. (2006) derived an analytical magnetic braking 
model assuming Ap/Bp stars for the formation of short period 
(P < 1 day) black hole low mass X-ray binaries and the AML 
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Figure 6. (a): The top panel shows the unfolded spectra along with the fit (thick line) for EPIC MOS1, MOS2 and PN. (b) The bottom panel displays the residuals. (c) 
The reference number is ObsID.0650020101 and observation date is 2010 December 11. 


is given by 


А 2 К 
Аль = GM) ABIIT i 


wind? (3) 
assuming M, = 0.52 Mo, R=0.88 Ro, B, (surface magnetic 
field) = 1000 G (upper limit) and Mwina = 107° Ms yr`!; Jii 
was estimated to be —2.54 x 10°° g cn? s ? and for Mind 
= 10-7? Ms уг! upper limit, Amb = —8.02 x 195 g cm? 52. 
RW Dor's spectral properties are similar to another source, VW 
Cep, which has been extensively studied both in optical and 


X-ray bands (Mitnyan et al. 2018). The assumption of magnetic 


10 


field B= 1000 G was based on observations of X-ray activity 
of VW Cep (Sanz-Forcada et al. 2007). Even a higher 
subsurface magnetic field of B — 20 kG is needed in VW Cep 
to explain the period modulation (Mitnyan et al. 2018). The 
angular momentum for a binary system can be written as 


q (GM?a) /?, 


= 4 
(1 + 4)? o 


о 
where q is mass ratio 0.63, М is the total mass of the system 
and a is the semimajor axis in cgs units resulting in 
Jo = 3.18 x 10?! in cgs units. Based on the period variation, 
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Figure 7. (а) The top panel shows the unfolded spectra along with the fit (thick line) for EPIC MOS1 and MOS2. (b): The bottom panel displays ће residuals. The 
reference number is ObsID.0723650201 and observation date is 2013 October 22. 
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we estimated P/P=9.81 x 10 ? yr ! and the AML rate is 
found to be J, = ZJ =3.31 x 10°° g cm? s ?, This value is a 


few times the value of Jm». This indicates that the stellar wind 
also plays a key role in removing the angular momentum from 
the system and helps in decreasing the orbital period of the 
system. The location of the spot over the stellar surface also 
affects the AML in the system. In a simulation, it has been 
found that higher latitude spots can trigger relatively more 
AML when compared to a spot located at lower latitude (Cohen 
et al. 2009). Moreover, it was also noticed that there was a 
feedback mechanism, i.e., AML affects the stellar dynamo as it 
controls the stellar magnetic field distribution. In our present 
study, we noted the spot located at a colatitude of 32? which is 
closer to the pole of the star and hence would be extracting 
relatively more AML. For such a study in contact binaries, it is 
necessary to know the role of spot location and AML. 

In the TRO model, as the material is transferring from 
primary to secondary, the decrease in the period should go 
along with an increase in the degree of contact, i.e., the CCE 
must be thick. However, we noted it is just 1096 which is 
marginal contact in nature. There could be two possible reasons 
for such a configuration. The first one is Liu et al. (2019) 
proposed a mechanism where the period of the contact binaries 
varies and hence the evolution of q, from higher to lower 
values, is driven by the oscillatory nature of fill-out factor f 
which depends on the thickness of the CCE. It was shown that 
the increase or decrease of f can decrease the period of the 
system and would be able to trigger a merger of the system. In 
the present study, f was close to 10%, i.e., CCE is not thick, and 
the mass transfer from the CCE to either of the stellar 
components has been ongoing as there is not sufficient material 
in the CCE. Perhaps, as the f continues to decrease, the material 
in the CCE would be transferred to the accepter star, i.e., the 
secondary component which is the second phase in the model 
of Liu et al. (2018). Another phenomenon could be mass loss 
from the CCE which can drive away some of the angular 
momentum and decrease the orbital period. 

It can be seen from O — C that, initially, the period was 
increasing and later it was driven to a period decreasing trend. 
The increasing trend can be explained by assuming a 
conservative mass transfer from the secondary to the primary 
component, however, there could be a mass loss from one or 
both of the components. The mass loss can be estimated from 
the following equation (Tout & Hall 1991) 

u -MMP 


2 P x 


assuming that the rate of mass transfer and loss is the same. We 
can equate M, from Equation (2) and the above results in the 
quadratic equation in terms of q, i.e., 34^ + 2q — 3 = 0, with 
one of the roots being q ~ 0.72. This indicates that at the high 
mass ratio, both mechanisms can drive the contact binary and 
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this also can explain the low degree of contact as some mass is 
lost from the system. 

Since RW Dor is displaying stellar activity in the form of the 
O'Connell effect, it is possible that the AML mechanism is 
occurring due to mass loss, which further helps in decreasing 
the orbital period. We assume that one of the components is 
driving the AML via mass loss, a situation which is described 
by Tout & Hall (1991) 

м = es ijs 6) 
2(M, + М) \ Ra) P 


1 


Here Ад is the Alfvén radius and d is the distance between the 
binary components. An upper limit on Alfvén radius can be 
estimated by assuming 2 = (d/RAY,, and conservative mass 
transfer and mass loss due to AML is the same, then it can be 
shown that q = M5/M, = 38 — 2/30 + 2 and for ~ 0.60, Ra 
= d/1.63 = 7.9 x 10'? cm which is similar to the size of the 
primary component radius. 

A detailed study of the O — C diagram was performed and 
light travel time effect (LITE) solutions with different epochs 
unveiled a third body with period Р; = 47.01 + 0.52 yr 
(Figure 4, top panel), which is similar to the results reported 
by Sarotsakulchai et al. (2019) but with zero eccentricity. In our 
analysis, the eccentricity was found to be ез = 0.21 which was 
obtained after performing a systematic search of the best e3 
based on the least residuals. Since a few of the times of minima 
were old and may not be accurate, we estimated the possibility 
of a third body's orbital period by removing the times of 
minima before 1980 and noted a period of P3 = 27.22 + 3.12 yr 
(see Figure 4, bottom panel). Such an orbital period of the third 
body was observed in VW Cep of about 31 yr (Kaszás et al. 
1998), which has very similar properties except for the mass 
ratio. Another source, BH Cas, which exhibits similar proper- 
ties to RW Dor and VW Cep also has a third body with 
P4 — 20.09 yr (Liu et al. 2019). For RW Dor, the mass of the 
third body was noted to be M3 = 0.175 Mo for i = 90°. We also 
looked for the third light in the TESS light curve but the data do 
not indicate any presence of it. Similar results were obtained in 
the previous study. Such a long period of third bodies is not 
rare and many of the contact binary systems have a third body 
with P4 of more than 20yr (Pribulla & Rucinski 2006; 
Rucinski & Pribulla 2007), and some of them have been 
optically confirmed. Future studies are needed to constrain the 
nature of third body solutions in RW Dor. 

Generally, such a periodic modulation is ruled out based on 
the Applegate mechanism (Applegate 1992) because it 
requires quadruple momentum AQ < 10°'-10°? (Lanza & 
Rodonó 2002). Lanza & Rodono (1999) showed that magnetic 
energy is proportional to the quadruple moment perturbations. 
However, a more recent theory by Lanza (2020) where the 
orbital angular momentum is turned into spin angular 
momentum notes that the AQ is larger by a factor of 
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Table 4 
Spectral Parameters Using a Bremsstrahlung Model and all the Error Bars at 90% Confidence Level 

Parameter Mos1(2009) Mos1(2010) Mos2(2010) PN(2010) Mos1(2013) Mos2(2013) 
kT (keV) 0.59 + 0.12 0.59 + 0.09 0.65 + 0.14 0.50 + 0.09 0.48 + 0.09 0.45 + 0.11 
Norm (x 10~*) 2.97 :E 1.12 3.51 + 1.06 3.49 + 0.92 2.21 + 0.48 3.9] + 1.15 3.27 + 1.09 
x7 /(dof) 45/15 36/18 29/15 59/41 27/18 19/11 

Table 5 

Spectral Parameters for bbody 

Parameter Mos1(2009) Mos 1(2010) Mos2(2010) PN(2010) Mos 1(2013) Mos2(2013) 
kT (keV) 0.18 + 0.01 0.16 + 0.02 0.17 + 0.01 0.18 + 0.01 0.17 + 0.01 0.17 + 0.02 
Norm, (x10~°) 3.64 + 0.49 3.42 + 0.65 3.65 + 0.56 2.80 + 0.29 3.90 + 0.65 3.22 + 0.60 
Flux (0.3-2.0 keV) 2.05 1.71 1.90 1.56 2.06 1.70 
x7 /(dof) 26/15 28/18 16/15 46/41 20/18 8/11 
Note. Fluxes are in the unit 107 ^ erg cm ? s~' and all the error bars are at the 90% confidence level. 

Table 6 

Spectral Parameters for Apec Model 

Parameter Mos1(2009) Mos 1(2010) Mos2(2010) PN(2010) Mos 1(2013) Mos2(2013) 
kT, (keV) 0.17 + 0.07 0.30 + 0.15 0.28 + 0.11 0.24 + 0.08 0.27 + 0.14 0.29 + 0.13 
Z 0.44 + 0.36 0.16 + 0.11 0.09 + 0.11 0.12 + 0.08 0.11 + 0.08 0.15 + 0.07 
Norm, (x10~*) 0.83 + 0.66 0.79 + 0.57 2.81 + 0.55 2.19 + 1.13 3.20 + 1.60 1.41 + 1.30 
kT> (keV) 0.87 + 0.11 0.92 + 0.23 0.82 + 0.43 0.93 + 0.15 1.03 + 0.32 1.08 + 0.44 
Norm,(x 107+) 0.96 + 1.10 3.24 + 1.65 2.34 + 1.15 1.66 + 0.65 1.84 + 0.85 1.93 + 0.81 
Flux (0.3-2.0 keV) 1.83 1.80 1.65 1.61 2.15 1.83 
Lx (10? erg s!) 3.33 3.28 3.00 2.93 3.92 3.33 
EM,(x 10°? ст °) 1.51 1.44 5.12 4.00 5.84 2.57 
ЕМ, (x10 cm?) 1.75 5.91 4.27 3.03 3.35 3.52 
xX (dof) 13/12 23/15 16/11 16/55 19/15 5/8 


2 1 


Note. Fluxes are in the unit 107"? ergcm ? s^ 


100-1000. Sarotsakulchai et al. (2019) explained the periodic 
variations using the presence of a third body and ruled out the 
magnetic activity. With the new theory put forward by Lanza 
(2020), now the magnetic activity cycle can be one of the 
viable mechanisms for the observed periodic modulation as 
observed in this source. 


6.3. Loop Size Emitting the Quiescent X-Ray Emission in 
RW Dor 


For a few contact binaries coronal emission temperatures are 
known. In Table 4, we show the X-ray spectral parameters of 
four contact binaries including RW Dor from the present work. 
Three contact binaries, viz. RW Dor, VW Cep and 44 Boo, 
have similar orbital periods and spectral types. Among these 
four systems, BH Cas has the highest X-ray luminosity, i.e., 
Ly = 92.5 x 10? erg 5! and others have similar Lx (Table 4). 
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and all the error bars are at the 90% confidence level. 


Based on present studies, it can be concluded that during the 
X-ray observations by XMM-Newton, RW Dor and other 
contact binaries (VW Cep; Gondoin 2004a, 44 Boo; 
Gondoin 2004b, and BH Cas; Liu et al. 2019) were in a 
quiescent state for some duration. The exact reason for coronal 
X-ray emission in a quiescent state is debatable in terms of the 
distribution of magnetic field structures but the presence of 
loops cannot be ruled out. Ventura et al. (1998) explained the 
quiescent emission using two distinct classes of loops 
associated with the corona in G-type stars. One of the coronal 
loops is cooler with a temperature Tmax ~ 1.5—5 x 10° К and 
the other опе is relatively hotter with a temperature 1-3 x 
10’ K. A detailed loop modeling is not possible for the present 
work but the observed temperature (see Table 7) from the 
X-ray spectra does indicate the presence of such loops 
primarily based on the coronal temperatures, viz. 0.27 and 
1.03 keV (year 2013) (3.1 and 12 MK). 
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Table 7 
Comparison of X-Ray Spectral Parameters of Four Contact Binaries along with 
Other Parameters 


Parameter RW Dor VW Cep BH Cas 44 Boo 
Mitnyan et al. 
Present (2018), Gon- Liu et al. Gondoin 
work doin (2004a) (2019) (2004b) 
kT, (keV) 0.28 0.26 0.13/ 0.25 
0.26" 
kT, (keV) 0.82 0.56 © 0.65 
Lx 10” erg 57! 3.3 6.28 92.5 2.24 
P (days) 0.285 0.278 0.405 0.267 
q (mass ratio) 0.61 0.30 0.47 0.49 
P (day уг!) 2.8 | 18.5 | 32.7 1 E 


Notes. Up and down arrow symbols signify increase and decrease in period 
respectively. 
* Estimated from blackbody and bremsstrahlung models. 


Assuming that the loop systems are not dynamic and each 
loop consists of similar nonvarying pressure p, temperature 
T (K) and cross-section, then a loop length scale (L) can be 
obtained (Mewe et al. 1982) using the relation T — 1400(pL)!/ ? 
(Rosner et al. 1978) 

Lio = 74 x F x T? x ЕМУ х (А/В), (7) 
where Г.о is the loop half length in units of ior? cm, 77 is the 
coronal region temperature in units of 10’ К, EMs> is the 
emission measure in units of 10°? cm ? and F is the filling 
factor. Actually, F is a difficult parameter to calculate as it 
depends on many parameters of the loop, for example, the 
area of the loop. In the case of RW Dor, it is not known, 
however for 44 Boo, F was reported to be 5096—7096 
(Gondoin 2004b). Adopting a value of F = 0.7, R= 0.88 Ro 
and the rest of the values of parameters from Table 7, we 
estimated L= 1 x 10'°cm for temperature 0.93 keV (PN 
2010). For other coronal temperatures, it varies from 0.1 to 
6 x 10? cm. These values can be considered as upper limits 
due to uncertainty in the value of F. It can be observed that 
the estimated loop size from X-ray spectra is almost 7-8 
times lower than that estimated from the Alfvén radius. 
There are many approximations in estimating the loop length 
as discussed with one of them being the constant loop cross- 
section. If one considers a non-uniform cross-section of the 
order of 10 (Schrijver 1987; with 1 being a uniform cross- 
section) then the loop size increases. The estimated loop 
sizes are close to the sizes seen in 44 Boo and VW Cep 
(O'Connell 1951; Gondoin 2004b). Perhaps the estimated 
smaller loop size could be due to the fact that during the 
X-ray observations, the source was in a quiescent state, 
hence the estimated loop size was smaller. Now during the 
TESS observations, it is probably higher due to enhanced 
activity. 
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7. Conclusions 


Optical and X-ray studies have been carried out for a 
marginal contact binary RW Dor using TESS and XMM- 
Newton telescopes. For the first time, we report the variable 
O’Connell effect in the light curve and the maximum 
asymmetry in the light curve can be explained by a stellar 
spot of size 15° on the primary component at HJD 
2458489.2509. Using new times of minima, we produce a 
robust estimate of the third body orbital period P3 ~ 47 yr with 
an eccentricity ез = 0.21. We show that non-conservative mass 
loss plays a key role in the evolution of the binary apart from 
conservative mass transfer from the primary to the secondary 
component. The X-ray luminosity did not vary significantly 
during three different occasions of observation, indicating a 
constant X-ray quiescent emission from the source. Assuming 
that the quiescent emission arises from an undisturbed loop, its 
size was estimated to be around 0.6-1.0 x 10' cm. Future 
observations in optical and X-ray energy bands would be 
helpful to understand the nature of coronal geometry in RW 
Dor and the associated activity. 
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